Audio-based caricature exaggeration

ABSTRACT

A method that uses at least one hardware processor for receiving a three-dimensional model of an object, receiving an audio sequence embodied as a digital file that comprises a musical composition, generating a video frame sequence, and synthesizing the audio sequence and the video frame sequence into an audiovisual clip. The three-dimensional model is embodied as a digital file that comprises a representation of the object. The generating step comprises computing a caricature of the object by applying a computerized caricaturization algorithm to the three-dimensional model. The computing has scaling gradient fields of surface coordinates of the three-dimensional model by a function of a Gaussian curvature of the surface, and finding a regular surface whose gradient fields fit the scaled gradient fields. The computing is with a different exaggeration factor for each of multiple ones of the video frames, and the different exaggeration factor is based on one or more parameters of the musical composition of the audio sequence.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims the benefit of U.S. Provisional Patent Application No. 61/976,510, entitled “Audio-Based Caricature Exaggeration”, filed Apr. 8, 2014, the contents of which are incorporated herein by reference in their entirety.

BACKGROUND

A caricature is often defined as an illustration of an object (e.g. a person, an animal or even an inanimate object) in which some features are exaggerated (also “enhanced” or “emphasized”) while others are de-emphasized (also “oversimplified”). Artistically, caricatures are many times targeted at creating a comic or grotesque effect.

Some of the earliest caricatures are found in the works of Leonardo da Vinci, who actively sought people with deformities to use as models. The point was to offer an impression of the original which was more striking than a traditional portrait.

Caricaturization, traditionally, has been a wholly-manual task performed by skilled artists or hobbyists. However, the past few decades have seen many attempts to computerize the process of converting an image or a model of a certain object into a caricature.

The foregoing examples of the related art and limitations related therewith are intended to be illustrative and not exclusive. Other limitations of the related art will become apparent to those of skill in the art upon a reading of the specification and a study of the figures.

SUMMARY

The following embodiments and aspects thereof are described and illustrated in conjunction with systems, tools and methods which are meant to be exemplary and illustrative, not limiting in scope.

One embodiment relates to a method comprising using at least one hardware processor for: receiving a three-dimensional model of an object; receiving an audio sequence; generating a video frame sequence, said generating comprising computing a caricature of the object, wherein (a) the computing is with a different exaggeration factor for each of multiple ones of the video frames, and (b) the different exaggeration factor is based on one or more parameters of the audio sequence; and synthesizing the audio sequence and the video frame sequence into an audiovisual clip.

Another embodiment provides a computer program product comprising a non-transitory computer-readable storage medium having program code embodied thereon, the program code executable by at least one hardware processor for: receiving a three-dimensional model of an object; receiving an audio sequence; generating a video frame sequence, said generating comprising computing a caricature of the object, wherein (a) the computing is with a different exaggeration factor for each of multiple ones of the video frames, and (b) the different exaggeration factor is based on one or more parameters of the audio sequence; and synthesizing the audio sequence and the video frame sequence into an audiovisual clip.

Yet a further embodiment provides a system comprising: (1) a non-transitory computer-readable storage medium having program code embodied thereon, the program code comprising instructions for: receiving a three-dimensional model of an object, receiving an audio sequence, generating a video frame sequence, wherein the generating comprises computing a caricature of the object, wherein (a) the computing is with a different exaggeration factor for each of multiple ones of the video frames, and (b) the different exaggeration factor is based on one or more parameters of the audio sequence, and synthesizing the audio sequence and the video frame sequence into an audiovisual clip; and (2) at least one hardware processor configured to execute the instructions.

In some embodiments, the method further comprises using the at least one hardware processor for determining the one or more parameters for each of multiple periods of the audio sequence.

In some embodiments, the one or more parameters are selected from the group consisting of: amplitude, frequency and tempo.

In some embodiments, the program code is further executable by said at least one hardware processor for determining the one or more parameters for each of multiple periods of the audio sequence; and the one or more parameters are selected from the group consisting of: amplitude, frequency and tempo.

In some embodiments, the generating further comprises altering a view angle of the caricature along the video frame sequence.

In some embodiments, the three-dimensional model comprises a polygon mesh representation of the object.

In some embodiments, the exaggeration factor is applied uniformly, to the entirety of the three-dimensional model.

In some embodiments, the exaggeration factor is applied non-uniformly, only to one or more portions of the three-dimensional model, which portions amount to less than the entirety of the three-dimensional model.

In some embodiments, the computing of the caricature of the object comprises: constructing a look-up table comprised of (a) different visualizations of the caricature, each computed with one of the different exaggeration factors, and (b) the exaggeration factor for each of the different visualizations; and using each caricature visualization from the look-up table when the exaggeration factor of that caricature visualization is determined to be suitable for the one or more parameters of the audio sequence.

In some embodiments, the computing of the caricature of the object comprises: scaling gradient fields of surface coordinates of the three-dimensional model by a function of a Gaussian curvature of the surface; and finding a regular surface whose gradient fields fit the scaled gradient fields.

In some embodiments, the computing of the caricature of the object further comprises amplifying the scaling according to local discrepancies between the object and a reference object.

In some embodiments, the computing of the caricature of the object further comprises amplifying the scaling according to local discrepancies between the object and a scaled down version of the object.

In addition to the exemplary aspects and embodiments described above, further aspects and embodiments will become apparent by reference to the figures and by study of the following detailed description.

BRIEF DESCRIPTION OF THE FIGURES

Exemplary embodiments are illustrated in referenced figures. Dimensions of components and features shown in the figures are generally chosen for convenience and clarity of presentation and are not necessarily shown to scale. The figures are listed below.

FIG. 1 shows a depiction of gradient fields of coordinate functions at each point on a surface residing on a tangent plane;

FIG. 2 shows the results of exemplary caricaturization of a facial 3D model, based on a reference (also “average”) model;

FIGS. 3A-D visualize four of the equations of the Implementation Considerations section below;

FIG. 4 shows the results of exemplary caricaturization of an armadillo 3D model, with a visualization of the approximated Gaussian curvature of the original mesh;

FIG. 5 shows the results of exemplary caricaturization of a facial 3D model, with a visualization of the approximated Gaussian curvature of the original mesh;

FIG. 6 shows the results of exemplary caricaturization of a textured face 3D model;

FIG. 7 shows the results of exemplary caricaturization of another facial 3D model;

FIG. 8 shows a histogram of relative geodesic distortion (in %) between caricatures of a Centaur;

FIG. 9 shows a histogram of relative geodesic distortion (in %) between caricatures of a cat.

FIG. 10 shows the results of exemplary caricaturization of a centaur 3D model, in two different poses;

FIG. 11 shows the results of exemplary caricaturization of a dog 3D model, in four different poses;

FIG. 12 shows the results of exemplary caricaturization of a gymnast 3D model, in four different poses;

FIG. 13 shows the results of exemplary caricaturization of a martial artist 3D model, in four different poses;

FIG. 14A shows the results of Multidimensional Scaling of geodesic distance matrices for uncaricaturized 3D shapes;

FIG. 14B shows the results of Multidimensional Scaling of geodesic distance matrices for caricaturized versions of the 3D shapes of FIG. 14A;

FIG. 15A shows the results of Multidimensional Scaling of geodesic distance matrices for additional uncaricaturized 3D shapes of; and

FIG. 15B shows the results of Multidimensional Scaling of geodesic distance matrices for caricaturized versions of the 3D shapes of FIG. 15A.

DETAILED DESCRIPTION

Disclosed herein is a method for audio-based caricature exaggeration. The method may result in an audiovisual clip having a video channel and an audio channel; the audio channel may include an audio sequence, whereas the video channel may include a video frame sequence depicting a caricature of a certain object, wherein the caricature visually changes over time based on one or more audio parameters of the audio sequence.

The visual change to the caricature along the audiovisual clip may be bidirectional or unidirectional. A bidirectional change may include alternating periods of increasing and decreasing exaggeration of the caricature along the audiovisual clip. Conversely, a unidirectional change may include either an increasing or a decreasing exaggeration of the caricature along the entire length of audiovisual clip.

Such temporal changes to the visualization of the caricature may include a change to an exaggeration of one or more features of the object, and/or a change to a de-emphasis of one or more features of the object. Merely as an illustrative example, the resulting audiovisual clip may include a caricature of a woman's face, where the size of her nose and ears changes with some correspondence to audible changes in the music of the clip.

Accordingly, the method for audio-based caricature exaggeration, which yields an audiovisual clip, may include some or all of the following steps:

In one step of the method, a three-dimensional (3D) model of an object may be received. Examples of such an object include a whole or a partial human body (e.g. a human head, a complete human figure, etc.), an animal figure, a figure of an inanimate object, etc. The 3D model may be received, for example, as a digital file including a 3D polygon (e.g. triangle) mesh representing the object, as known in the art.

In another step of the method, an audio sequence may be received, in the form of a digital audio file. The audio sequence may include, for example, a musical composition, a song, speech, or any combination thereof. The audio sequence may serve as an audio channel of the audiovisual clip.

A further step of the method may include generating a video frame sequence, to serve as a video channel of the audiovisual clip. This step may include multiple discrete computations performed on the 3D model, each yielding a caricature visualization having a different level of exaggeration. Each of these computations may include applying a computerized caricaturization algorithm to the 3D model. An exemplary suitable algorithm, referred to as “Intrinsic Surface Caricaturization”, is disclosed below. However, those of skill in the art will recognize that other computerized caricaturization algorithms may be suitable for yielding caricature visualizations with different levels of exaggeration. Generally, any computerized caricaturization algorithm which allows adjusting the level of exaggeration of the resulting caricature may be suitable for the purpose of generating the present video frame sequence.

In the Intrinsic Surface Caricaturization algorithm, adjusting the level of exaggeration may be achieved by adjusting an exaggeration factor, termed γ (Gamma). Namely, in each of the multiple discrete computations performed on the 3D model, a different γ value (γε[0,1]) may be used.

The γ value may be adjusted either for the entirety of the 3D model, or for one or more specific portions thereof, which amount to less than the entirety of the model. In the latter case, the object may be exaggerated non-uniformly; one or more specific portions of the object may be exaggerated while others remain static. This may be achieved by setting the exaggeration factor non-uniformly over the 3D model of the object; a portion which is desired to remain static may have an exaggeration factor of zero or close to zero, while the rest of the 3D model may have an exaggeration factor relatively distant from zero, e.g. 0.25, 0.5, 0.75, 1, etc. With regard to the Intrinsic Surface Caricaturization algorithm, those of skill in the art will recognize that γ values may be differently set for different portions of the 3D model.

Optionally, the user may manually specify the portion(s) of the 3D model to be exaggerated, such as by marking, using a graphical user interface (GUI), one or more points, triangles or the like in the 3D model.

Practically, the generation of the video frame sequence may include a preliminary process of constructing a look-up table comprised of different caricature visualizations and the γ value for each. The look-up table may be predetermined to include a certain number of caricatures, e.g. 10 caricatures, 50 caricatures, 100 caricatures, etc. The number of caricatures in the look-up table may, naturally, imply the γ value increments of the multiple computations. However, these increments are not necessarily equal, and it is possible, if desired, to produce caricatures with a non-liner γ value incrementation.

After the look-up table is constructed, the caricatures it contains may be conveniently used for generating the frames of the video frame sequence, by using a certain caricature visualization from the look-up table when its γ value is determined to be suitable for one or more audio parameters of the audio sequence, at a certain point in time or during a certain time period.

As an alternative to constructing a look-up table, a computation with a desired γ value may be performed ad-hoc, for each video frame, based on the one or more audio parameters corresponding to that frame. If this approach is used, a caricature with a certain γ value may be re-used for a later frame if needed, without repeating the same computation for that later frame.

In the generation of the video frame sequence, the level of caricature exaggeration may be synchronized with the audio sequence in the following manner: For each period of time in the audio sequence (e.g. a few milliseconds, a few dozen milliseconds, a few hundred milliseconds, or any different period of time defined by a user), one or more audio parameters of that period in the audio sequence may be determined, by analyzing the audio sequence. Such audio parameters may include, for example:

A. Amplitude (also “volume”) of the audio sequence during the pertinent period of time. This may be, for example, a dominant, median or average amplitude over the pertinent period of time. The amplitude may relate to the entire frequency spectrum present during the pertinent period of time, or to a specific band of frequencies.

B. Sound frequency (also “pitch”) of the audio sequence during the pertinent period of time. This may be, for example, a dominant, median or average frequency over the pertinent period of time.

Those of skill in the art will recognize that any other human-audible feature of the audio sequence may serve as such a parameter, such that the caricature changes upon the human hearing such feature.

In addition or as an alternative to determining the one or more audio parameters for every period of time in the audio sequence, a global audio parameter pertaining to the entirety of the audio sequence (or to an extended section thereof, e.g. 10 seconds or more) may be determined, and be used for the synchronization. An example of such global audio parameter is tempo; the tempo, denoted, for example, in beats per minute (BPM), may be used for determining a frequency of changing the visualization of the caricature. The frequency of chancing the visualization of the caricature may be determined as an integer multiple of the BPM. For example, if the BPM is 100, the visualization of the caricature may be changed 100 times per minute, 200 times per minute, 300 times per minute, and so on and so forth.

After determining the one or more audio parameters for every period of time in the audio sequence, one or more video frames may be generated, to be associated with that period of time. For example, if it is desired to produce an audiovisual clip having a frame rate of 25 frames per second (FPS), which translates into 40 milliseconds per frame; and the length of each of the aforementioned periods of time is 80 milliseconds, then two frames have to be generated for cach such period of timc.

As noted above, the one or more video frames which are generated per period of time in the audio sequence may include a caricature visualization having a level of exaggeration which corresponds to a level of the one or more audio parameters. To that end, the audio sequence may be pre-analyzed, to determine how the level of the one or more audio parameters changes during the audio sequence. These changes may then be encoded on a scale which correlates to the scale of the exaggeration factor, for example between 0 and 1. This may make it convenient to select a level of exaggeration which corresponds to the level of the one or more audio parameters.

It should be noted that the computation of a caricature visualization need not be necessarily carried out for all frames in the video channel. It is possible to compute a caricature visualization only once every few frames, and generate the frames in between using morphing techniques known in the field of digital video.

It should also be noted that, in order for the resulting video channel to have a smooth, appealing appearance to a viewer, the computation of the caricature visualizations of different levels may be performed not based on the audio sequence itself, but rather on a smoother version thereof. That is, the audio sequence may be averaged, for example by using a low pass filter, such that the changc in the one or morc audio paramctcrs along the averaged version of the audio sequence becomes morc subtle (i.e. of a slower pace). This will cause the resulting video channel to include caricature visualizations whose temporal changing is smoother and less bothering.

In an optional step of the method, a view angle of the caricature along the video frame sequence may also be changed, either based on the one or more audio parameters or based on some user preference unrelated to the one or more audio parameters. To achieve this, the caricaturization algorithm may be separately applied to different postures of the 3D model of the object. The result may be, for example, a video channel in which the caricature both changes its exaggeration level and is shown from different angles (e.g. if the caricature is of a human head, the head may also rotate during playing of the video channel).

In a further step of the method, the audio sequence and the video frame sequence may be synthesized into an audiovisual clip. Namely, the audio sequence and the video frame sequence may be joined together, into a single file, as known in the art.

In an optional step of the method, the audiovisual clip may undergo post-processing, for example in order to compress its audio sequence and/or it video frame sequence. This optional step may employ one or more of audio and/or video encoding algorithms known in the art.

The audiovisual clip produced according to the present method may show smooth, appealing transformations between such caricature visualizations having different exaggeration factors (and/or between an original and caricatures of the original)—while the audio sequence is playing in the background. One possible use of such audiovisual clips is for entertainment purposes. Other possible uses will become apparent to those of skill in the art.

Intrinsic Surface Caricaturization

1. Surface Enhancement

1.1. Notations

Consider the boundary of a 3D shape as a surface or a two dimensional manifold M⊂

³ given in a parametric form as S: Ω

²→M. Throughout this exposition, the surface M will be described by its coordinates S(u, v)={x(u, v), y(u, v),z(u, v)}, where x, y and z are the coordinate-functions of the manifold.

The first partial derivatives of the parametric form, S_(u) and S_(v), are two vectors that are tangent to the surface and linearly independent. Thus, one can set the basis for the tangent plane T_(p)M, about a point pεS on the surface, as {S_(u),S_(v)}.

The Jacobian of the parametrization J is a 3×2 matrix whose columns are S_(u) and S_(v), respectively. We define the matrix G=J^(T) J as the regular metric of the surface. Clearly, G is invertible, positive definite, and can be written explicitly in terms of S_(u) and S_(v) as

$\begin{matrix} {G = {\begin{pmatrix} {\langle{S_{u},S_{u}}\rangle} & {\langle{S_{u},S_{v}}\rangle} \\ {\langle{S_{v},S_{u}}\rangle} & {\langle{S_{v},S_{v}}\rangle} \end{pmatrix}.}} & (1) \end{matrix}$

Since the basis {S_(u),S_(v)} is not necessarily orthogonal, the inner product between two vectors a, b εT_(p)M represented in this basis is given by

a,b

_(G) =a ^(T) Gb.  (2)

Next, consider two surfaces M₁ and M₂ and suppose we are given a parametrization for the first surface S₁(u,v). Suppose we are given a bijective map between the surfaces φ:M₁→M₂. We can then build a parametrization for the second surface M₂ by composing φ on S₁. That is, a parametrization for the second surface is given by

S ₂(u,v)=φ∘S ₁(u,v).  (3)

Given the parametrizations for both surfaces, we can calculate the metric matrices, G₁ and G₂ at each point (u, v) in the parametrization domain, see Eq. (1). If at each point (u, v) in the parametrization domain, G₁ (u, v)=G₂ (u, v), then φ is said to be an isometry. If there exists an isometry between two surfaces M₁ and M₂ they are said to be isometric.

A curve on one surface has the same length as that of its isometrically mapped version. An area of a patch on one surface is mapped by φ to a patch that has the same area. In fact, any geometric property that can be measured in terms of the metric G and its derivatives and integration with respect to u and v are preserved by the isometry φ. We call such properties isometric invariant or intrinsic.

For each function ƒ: S→

defined on the surface there exists a corresponding function in the parameter domain {tilde over (ƒ)}: Ω⊂

²→

, such that

{tilde over (ƒ)}(u,v)=ƒ∘S(u,v).  (4)

Similarly, we define a tangent vector field on the surface, {right arrow over (F)}: M→T_(p)M, as an assignment of a tangent vector to each point on the surface. Since these fields are in the tangent plane, one can represent them in terms of the basis {S_(u),S_(v)}. For each tangent vector field {right arrow over (F)} on the surface, there exists a corresponding vector field with identical coordinates in the parametrization domain {tilde over ({right arrow over (F)})} represented by the natural basis of

².

As defined in M. P. Do Carmo, “Dierential geometry of curves and surfaces”, Vol. 2, Prentice-hall Englewood Cliffs, 1976, the gradient of a function on the surface, which we term here as gradient with respect to the metric G, is a differential operator that satisfies

$\begin{matrix} {{{\langle{{\nabla_{G}{f(p)}},\overset{\rightarrow}{w}}\rangle}_{G} = {{d_{p}f} = {\lim\limits_{h\rightarrow 0}\frac{{f\left( {p + {h\overset{\rightarrow}{w}}} \right)} - {f(p)}}{h}}}},} & (5) \end{matrix}$

for each {right arrow over (w)}εT_(p)M. It follows that

ƒ(p+h{right arrow over (w)})=ƒ(p)+h

∇ _(G)ƒ(p),{right arrow over (w)}

_(G) +O(h ²).  (6)

For the corresponding function in the parametrization domain {tilde over (ƒ)}(u,v),

ƒ(u ₀ +hw ₁ v ₀ +hw ₂)=ƒ(u ₀ ,v ₀)+h(∇_(u,v){tilde over (ƒ)}(u ₀ ,v ₀),{right arrow over (w)}

+O(h ²),  (7)

where p=S(u₀,v₀), and {right arrow over (w)}=(w₁, w₂). Since the left hand sides of (6) and (7) are equal, it follows that

∇_(G)ƒ(p),{right arrow over (w)}

_(G)=

∇_(u,v){tilde over (ƒ)}(u ₀ ,v ₀),{right arrow over (w)}

.  (8)

And, therefore,

{right arrow over (w)} ^(T) G∇ _(G)ƒ(p)={right arrow over (w)} ^(T)∇_(u,v){tilde over (ƒ)}(u ₀ ,v ₀),  (9)

for each {right arrow over (w)}εT_(p)M. This provides a formula for the gradient in terms of the gradient of its corresponding function in the parametrization domain {tilde over (ƒ)} and the parametrization S as

∇_(G)ƒ(p)=G ⁻¹∇_(u,v){tilde over (ƒ)}(u ₀ ,v ₀).  (10)

Treating the space of functions and the space of fields on the surface as an Hilbert space, one can define inner product between two functions, ƒ and g, as

ƒ,g

_(L) ₂ =∫_(Ω)ƒ(x)g(x)da  (11)

and the inner product between to vector fields, {right arrow over (F)} and {right arrow over (H)}, as

{right arrow over (F)},{right arrow over (H)}

_(F) ₂ =∫_(Ω) {right arrow over (F)}/(x)^(T) G{right arrow over (H)}(x)da.  (12)

Here, da=√{square root over (det G)}dudv.

The divergence is a differential operator that produces a function from a tangent field. It is defined as minus the adjoint operator of the gradient, with respect to the above inner products. That is, the divergence of fields on a manifold, denoted as ∇_(G′), satisfies

∇_(G) ·{right arrow over (F)},g

_(L) ₂ =−

{right arrow over (F)},∇ _(G) g

_(F) ₂ ,  (13)

for any function g and a field {right arrow over (F)} on the surface. By integration by parts and assuming that g vanishes on the boundaries, one can show that the divergence with respect to the metric of a tangent field can be calculated by

$\begin{matrix} {{{\nabla_{G}{\cdot \overset{\rightarrow}{F}}} = {\frac{1}{\sqrt{\det \mspace{14mu} G}}{\nabla_{u,v}{\cdot \left( {\sqrt{\det \mspace{14mu} G}\overset{\rightarrow}{\overset{\sim}{F}}} \right)}}}},} & (14) \end{matrix}$

where {tilde over ({right arrow over (F)})} is the corresponding field of {right arrow over (F)} in the parameterization domain.

The laplacian is a differential operator that maps functions to functions. The laplacian of functions on the surface is called the Laplace-Beltrami operator and is denoted here as Δ_(G). It is defined as the divergence of the gradient. Thereby, it can be written as

$\begin{matrix} {{\nabla_{G}{\cdot f}} = {\frac{1}{\sqrt{\det \mspace{14mu} G}}{\nabla_{u,v}{\cdot \left( {\sqrt{\det \mspace{14mu} G}G^{- 1}{\nabla_{u,v}\overset{\sim}{f}}} \right)}}}} & (15) \end{matrix}$

1.2 Gradient-Based Deformation

The gradient fields of these functions with respect to a given metric, ∇_(G)x, ∇_(G)y and ∇_(G)z, define three vector fields that reside in the tangent plane T_(p)M, defined by S_(u) and S_(v), about each point s on the surface pεM. See FIG. 1, which shows the gradient fields of the coordinate functions at each point on the surface residing on the tangent plane.

The gradient-based deformation method, proposed in Y. Yu, K. Zhou, D. Xu, X. Shi, H. Bao, B. Guo, H.-Y. Shum, “Mesh editing with poisson-based gradient field manipulation”, ACM SIGGRAPH 2004 vol. 23 (3), 644-651, first transforms each field locally and arbitrarily, ∇_(G)x(p)→T_(p)(∇_(G)x(p)). In order to derive the deformed surface, {tilde over (S)}={{tilde over (x)}, {tilde over (y)}, {tilde over (z)}}, whose coordinate-functions gradient fields mostly agree with the manipulated ones in a least-squares sense, the method minimizes the following term for each coordinate function separately

∫_(S)|∇_(G) {tilde over (x)}(p)−T _(p)(∇_(G) x(p))∥² da(p),  (16)

with similar functionals for y and z. Here, da(p) is an area element about the point pεS. The Euler-Lagrange equation of the above functional is the Poisson equation

Δ_(G) {tilde over (x)}=∇ _(G)·(T _(p)(∇_(G) x)),

or in surface notation,

Δ_(G) {tilde over (S)}=∇ _(G)·(T _(p)(∇_(G) S)).  (17)

This simple method gave birth to a rich family of surface deformation techniques. For example, restricting the deformations to isometrics, the set of local transformations {T_(p)} is limited to rotation matrices. In the present exaggeration method, we restrict ourselves to simple scalar scaling.

Other linear mesh deformation techniques include thin-shell (M. Botsch, L. Kobbelt, “An intuitive framework for real-time freeform modeling”, ACM Transactions on Graphics (TOG) 23 (3), 2004, 630-634) and Laplacian-based (O. Sorkine, D. Cohen-Or, Y. Lipman, M. Alexa, C. Rossi, H.-P. Seidel, “Laplacian surface editing”, 2004, 175-184), to name just a few. M. Botsch, O. Sorkine, “On linear variational surface deformation methods”, Visualization and Computer Graphics, IEEE Transactions on, 14 (1), 2008, 213-230, reviews some of these techniques and compares them. For a broader overview of these and other related deformation methods, see M. Botsch, L. Kobbelt, M. Pauly, P. Allicz, B. Levy., “Polygon mesh processing”, A K Peters/CRC Press, 2010.

Next, we define a model-based framework for relative exaggeration. Let S be the surface of a given object. In addition, assume we are given a bijective map φ: S→Ŝ, between S and the surface, Ŝ, of a given representative average object. Such maps can be obtained by surface registration algorithms like non-rigid ICP with a sparse correspondence (see T. Weise, H. Li, L. Van Gool, M. Pauly, “Face/off: Live facial puppetry”, in Proceedings of the 2009 ACM SIGGRAPH/Eurographics Symposium on Computer animation, ACM, 2009, 7-16).

In accordance with present embodiments, there may be two requirements from a caricaturization process: First, it needs to exaggerate the discrepancies from a given model. Second, it needs to be invariant to pose or independent of the embedding. To satisfy the latter, we scale the given surface gradient fields by some function of the absolute Gaussian curvature |K| of the surface S. The choice of the Gaussian curvature, the simplest isometry invariant differential quantity, as an amplification factor makes the present method robust to isometrics of the surface.

In order to exaggerate the discrepancies from the average model, we choose the following function as a local scaling factor

$\begin{matrix} {{{f(p)} = {{K(p)}}^{\beta \mspace{11mu} \log \frac{{da}{(p)}}{d{\hat{a}{({\varphi {(p)}})}}}}},} & (18) \end{matrix}$

where β is some positive constant, da(p) is an area element about pεS and dâ(φ(p)) is the corresponding area element about φ(p)εŜ— the average representative surface. Note that given φ, the corresponding area elements can be easily computed by defining a corresponding parametrization.

In warped surface elements where there is no local scaling, that is da(p)=da(φ(p)), no exaggeration is executed. Regions that were shrunk with respect to the reference model will shrink even more, while those that were stretched will stretch more in the present method. In the case where no prior average model exists, we could assume, for example, that log

$\frac{{da}(p)}{d{\hat{a}\left( {\varphi (p)} \right)}} = ɛ$

and thereby scale the gradient fields by |K|^(γ), where γ=εβ is some non-negative constant.

FIG. 2 shows, from left to right: a reference (also “average”) 3D model of a face; a 3D model of an actual face, which naturally differs from the reference; a caricature of the actual face, with an exaggeration factor β=0.3; and a caricature of the actual face, with an exaggeration factor β=0.5.

In order to find a regular surface whose gradients with respect to the parametrization best fits the scaled ones, we minimize the squared difference between the tangent fields and end up with the following Poisson equation

$\begin{matrix} {{\Delta_{G}\overset{\sim}{S}} = {{\nabla_{G}{\cdot \left( {{f(p)}{\nabla_{G}S}} \right)}} = {\nabla_{G}{\cdot {\left( {{{K(p)}}^{\beta \mspace{11mu} \log \frac{{da}{(p)}}{d{\hat{a}{({\varphi {(p)}})}}}}{\nabla S}} \right).}}}}} & (19) \end{matrix}$

Here, {tilde over (S)} is the exaggerated surface we are looking for, and the differential operators are taken with respect to the metric of the given surface S. When there is no reference surface Ŝ, the enhancement simplifies to the solution of the following Poisson equation

Δ_(G) {tilde over (S)}=∇ _(G)·(|K| ^(γ)∇_(G) S).  (20)

2. Isometry Preservation

We define two isometric surfaces as parametrized surfaces between which an isometric map exists. By construction, the present surface enhancement model is robust to isometric transformations. That is, given two isometric surfaces S₁S₂, if we can find corresponding surfaces {tilde over (S)}₁ and {tilde over (S)}₇ that satisfy Equation (20), we also have that {tilde over (S)}₁

{tilde over (S)}₂.

Thorem 2.1.

If the surfaces S₁ and S₂ are isometric, and there exist two surfaces {tilde over (S)}₁ and {tilde over (S)}₂ such that

∇_(G) ₁ {tilde over (S)} ₁ =|K ₁|∇_(G) ₁ S ₁,  (21)

and

∇_(G) ₂ {tilde over (S)} ₂ =|K ₂|∇_(G) ₂ S ₂,  (22)

then, {tilde over (S)}₁ and {tilde over (S)}₂ are isometric.

Proof.

By definition, the metric of {tilde over (S)}₁ can be written in terms of the Jacobian of the parametrization as

{tilde over (G)} ₁ ={tilde over (J)} ₁ ^(T) {tilde over (J)} ₁,  (23)

where {tilde over (J)}₁ is a matrix whose columns are {tilde over (S)}_(1u) and {tilde over (S)}_(1v).

By writing the gradient of the coordinate functions of {tilde over (S)}₁ as functions on the surface S₁ explicitly, we get

∇_(G) ₁ {tilde over (S)} ₁ =G ₁ ⁻¹ ∇{tilde over (S)} ₁ =G ₁ ⁻¹ {tilde over (J)} ₁.  (24)

Therefore,

{tilde over (G)} ₁=∇_(G) ₁ {tilde over (S)} ₁ ^(T) G ₁ G ₁∇_(G) ₁ {tilde over (S)} ₁.  (25)

By plugging relation (21) into the above equation, we get

{tilde over (G)} ₁ =|K ₁|²∇_(G) ₁ S ₁ ^(T) G ₁ G ₁∇_(G) ₁ S ₁ =|K ₁|² ∇S ₁ ^(T) ∇S ₁ =|K ₁|² G ₁.  (26)

By repeating the same derivation for {tilde over (S)}₂, we similarly obtain

{tilde over (G)} ₂ =|K ₂|² G ₂.  (27)

Since we assumed that the surfaces S₁ and S₂ are isometric, the metric and the Gaussian curvature are equal for corresponding points on the these surfaces. Hence,

|K ₁|² G ₁ =|K ₂|² G ₂.  (28)

We can conclude that

{tilde over (G)} ₁ ={tilde over (G)} ₂  (29)

which proves that indeed the surfaces {tilde over (S)}₁ and {tilde over (S)}₂ are isometric.

Since we minimize the difference between the fields on isometric surfaces, and the divergence is intrinsic, the exaggeration is robust to isometrics in a least squares sense. In practice, surfaces which are isometric to one another are rigidly equivalent in the discrete case. However, experimentally, for almost isometric shapes, the isometric deviation of the resulting caricatures, i.e. the Gromov-Hausdorff distance, is negligible.

3. Implementation Considerations

One possible implementation of the present framework is when approximate continuous surfaces by triangulated meshes. In what follows, we briefly review discrete approximations of a continuous function defined on the surface, its gradient, the surface Gaussian curvature, a discrete divergence operator (Y. Tong, S. Lombeyda, A. N. Hirani, M. Desbrun, “Discrete multiscale vector field decomposition”, Proceedings of ACM SIGGRAPH 2003, vol. 22 (3), 445-452), and a discrete Laplace-Beltrami operator (M. Meyer, M. Desbrun, P. Schroder, A. H. Barr, “Discrete differential geometry operators for triangulated 2-manifolds”, in Visualization and mathematics III, Springer, 2003, 35-57). We conclude by combining these elements together for solving the system of equations (20) with some additional constraints that would fix, for example, some surface points as boundary conditions.

For the discretization we use the finite element method (FEM), in which every function ƒ defined on the surface is assumed to be piecewise linear. As such, it can be represented with the piecewise linear basis functions {B_(i)}_(i=1) ^(n) as

${{f(x)} = {\sum\limits_{i = 1}^{n}{f_{i}{B_{i}(x)}}}},$

where n is the number of vertices, ƒ_(i) is the value of the function ƒ at the vertex v_(i), and B_(i)(x) is a linear B-spline function that is equal to 1 at the vertex v_(i) and 0 at the rest of the vertices. Within each triangle T=(v_(i),v_(j),v_(k)), ƒ is linearly interpolated only by ƒ_(i), ƒ_(j), and ƒ_(k). It can be shown that the gradient of ƒ within each triangle is constant (Meyer et al. 2003), and given by

${{\nabla{f(T)}} = {{\left( {f_{j} - f_{i}} \right)\frac{\left( {s_{i} - s_{k}} \right)^{\bot}}{2A_{T}}} + {\left( {f_{k} - f_{i}} \right)\frac{\left( {s_{j} - s_{i}} \right)^{\bot}}{2A_{T}}}}},$

illustrated in FIG. 3A, where s_(i), s_(j), and s_(k) are the coordinates of the vertices v_(i), v_(j), and v_(k), respectively, ⊥ denotes the counterclockwise 90° rotation operator in the triangle plane, and A_(T) is the area of the triangle T.

There are different computational schemes for locally estimating the Gaussian curvature of a triangulated mesh. See T. Surazhsky, E. Magid, O. Soldea, G. Elber, E. Rivlin, “A comparison of gaussian and mean curvatures estimation methods on triangular meshes”, Robotics and Automation, 2003. Proceedings. ICRA '03. IEEE International Conference on, (1), 2003, 1021-1026, for a review and evaluation. Introduced in Meyer et al. 2003 is the so called Gauss-Bonnet scheme, given by

${{K\left( v_{i} \right)} = {\frac{1}{A_{i}}\left( {{2\pi} - {\sum\limits_{j:{v_{j} \in_{1}{(v_{i})}}}\theta_{j}^{i}}} \right)}},$

illustrated in FIG. 3B, where A, is the Voronoi area about the vertex v_(i) and θ_(j) ^(i)s are the incident angles to the vertex v_(i). Another scheme, presented in D. Cohen-Steiner, J.-M. Morvan, “Restricted delaunay triangulations and normal cycle”, SCG '03 Proceedings of the nineteenth annual symposium on Computational geometry, 312-321, approximates the second fundamental form at each vertex v_(i). The eigenvalues of the second fundamental form are the principal curvatures, k₁ and k₂, by which we derive the Gaussian curvature to be K(v_(i))=k₁(v_(i))·k₂ (v_(i)). An implementation of this scheme in Matlab is given in Gabriel Peyre's Toolbox Graph. We tested both schemes in our experiments.

A discrete formulation of the divergence operator can be derived from the divergence theorem

∫∫_(A)

·{right arrow over (F)}da=

_(∂A) {right arrow over (F)}·{circumflex over (n)}dl,

where {right arrow over (F)} is a given vector field, and {circumflex over (n)} is an outward pointing unit vector normal to the boundary. We assume that {right arrow over (F)} is constant within each triangle of the mesh, and approximate its divergence at a vertex v_(i) by averaging the divergence over a 1-ring neighborhood,

${\begin{matrix} {{{{Div}\left( \overset{\rightarrow}{F} \right)}\left( v_{i} \right)} \cong \frac{\int{\int_{\;_{1}{(v_{i})}}{{\nabla{\cdot \overset{\rightarrow}{F}}}{da}}}}{\int{\int_{\;_{1}{(v_{i})}}{da}}}} \\ {= \frac{\sum\limits_{T_{j} \in_{1}{(v_{i})}}{\int{\int_{A{(T_{j})}}{{\nabla{\cdot \overset{\rightarrow}{F}}}{da}}}}}{\sum\limits_{T_{j} \in_{1}{(v_{i})}}{A\left( T_{j} \right)}}} \\ {= \frac{\sum\limits_{T_{j} \in_{1}{(v_{i})}}{\oint_{\partial{A{(T_{j})}}}{{\overset{\rightarrow}{F} \cdot \hat{n}}{dl}}}}{\sum\limits_{T_{j} \in_{1}{(v_{i})}}{A\left( T_{j} \right)}}} \\ {= \frac{\sum\limits_{T_{j} \in_{1}{(v_{i})}}{\sum\limits_{e_{k} \in T_{j}}{{F\left( T_{j} \right)} \cdot {e_{k}\left( T_{j} \right)}^{\bot}}}}{\sum\limits_{T_{j} \in_{1}{(v_{i})}}{A\left( T_{j} \right)}}} \end{matrix},}\quad$

illustrated in FIG. 3C, where _(i)(v_(i)) is the 1-ring neighborhood, comprised of the triangles incident to v_(i), A(T_(j)) is the area of the triangle T_(j), and e_(k) (T_(j))^(⊥) is an outwardly directed vector normal to the k^(th) edge of the triangle T_(j) with similar length to the edge e_(k)(T_(j)).

For discretizing the Laplace-Beltrami operator, we use the cotangent weight scheme for triangulated meshes (Meyer et al. 2003). We plug the values of a given function at each vertex, ƒ_(i)=ƒ(v_(i)), to a vector, and multiply it by L=A⁻¹W, where A is the diagonal matrix of the Voronoi cell areas about each vertex, as introduced in U. Pinkall, K. Polthier, “Computing discrete minimal surfaces and their conjugates”, Experimental Mathematics 2 (1), 1993, 1536. W is the classic cotangent weight matrix

$W_{ij} = \left\{ {\begin{matrix} {\sum_{j:{v_{j} \in {_{1}{(v_{i})}}}}\left( {{\cot \mspace{14mu} \alpha_{ij}} + {\cot \mspace{14mu} \beta_{ij}}} \right)} & {{{if}\mspace{14mu} i} = j} \\ {{{- \cot}\mspace{14mu} \alpha_{ij}} - {\cot \mspace{14mu} \beta_{ij}}} & {{{{if}\mspace{14mu} i} \neq j},{v_{j} \in {_{1}\left( v_{i} \right)}}} \\ 0 & {otherwise} \end{matrix},} \right.$

illustrated in FIG. 3D, where α_(ij) and β_(ij) are the angles opposite to the edge (v_(i), v_(j)).

Boundary Conditions:

There are various ways to deal with boundaries in the present framework. Here, we describe three possible alternatives. Other options could also be applied depending on the application. Next, we describe how to incorporate either soft or hard constraints along the boundary of the shape.

For hard constraints, we formulate the linear system as the constrained least-squares optimization problem

$\begin{matrix} \begin{matrix} \min\limits_{\overset{\sim}{x}} & {{{L\overset{\sim}{x}} - b}}_{A}^{2} \\ {s.t.} & {{x_{j} = x^{*}},{v_{j} \in {\partial S}},} \end{matrix} & (30) \end{matrix}$

where b=Div(|K|^(γ)Grad(x)). Note that the norm is weighted by the area of the Voronoi cells about any vertex. Therefore, the norm can be written as ∥F|_(A) ²=(F^(T)AF). The linear equality constraints can be written as B{tilde over (x)}=x*, where B is a sampling matrix, in which each row is zero, besides those that correspond to the constrained vertices. x* is a vector of the final positions of the constrained vertices. Using the Lagrange multipliers approach, the solution to the above problem is

$\begin{matrix} {{\begin{pmatrix} L & B \\ B^{T} & 0 \end{pmatrix}\begin{pmatrix} \overset{\sim}{x} \\ \lambda_{x} \end{pmatrix}} = {\begin{pmatrix} {{Div}\left( {{K}^{\gamma}{{Grad}(x)}} \right)} \\ x^{*} \end{pmatrix}.}} & (31) \end{matrix}$

Similar systems are solved for {tilde over (y)} and {tilde over (z)}.

For soft constraints, one can add a controllable quadratic term for manipulating the boundary vertices and minimize the following unconstrained optimization problem

$\begin{matrix} {{\min\limits_{\overset{\sim}{x}}{{{L\overset{\sim}{x}} - b}}_{A}^{2}} + {\mu {{{{B\overset{\sim}{x}} - x^{*}}}_{A}^{2}.}}} & (32) \end{matrix}$

The solution to Problem (32) is

{tilde over (x)}=(LAL+μB ^(T) AB)⁻¹(ALb+μAB ^(T) x*).  (33)

The third type of boundary constraints are set by fixing γ=0 on selected triangles. This technique enables keeping some regions fixed throughout the exaggeration.

3. Experimental Results

We evaluate the present caricature generation method by conducting and discussing several example use cases. First, we demonstrate its applicability by caricaturizing a number of synthetic surfaces and captured geometries. Then, we discuss the performance analysis of a real-time implementation of the present algorithm. Next, we examine the robustness of the method to almost isometric deformations. We conclude the section by proposing two possible applications.

3.1 Qualitative Evaluation

We demonstrated the power of the self-caricaturization method by applying it to triangulated meshes. The shade of the meshes in FIGS. 4 and 5 represents the approximated Gaussian curvature at each point.

FIG. 4 shows, from left to right: an original 3D model of a armadillo; a caricature with an exaggeration factor γ=0.5; and a caricature with an exaggeration factor γ=1. The approximated Gaussian curvature of the original mesh is indicated by the darkness of the armadillo, where larger values are highlighted using a darker shade, visible especially at the armadillo's fingertips, toe tips, ear tips, tail tip, nose tip, etc., and smaller values are shown as lighter shades.

FIG. 5 shows, from left to right: an original facial 3D model; a caricature with an exaggeration factor γ=0.5; and a caricature with an exaggeration factor γ=0.75. The approximated Gaussian curvature of the original mesh indicated by the darkness of the face, where larger values are highlighted using a darker shade, visible especially at the bottom part of the nose, the lips and the eyes, and smaller values are shown as lighter shades. Note that the darker shade of the face's boundaries does not imply that they have any substantial Gaussian curvature. In fact, areas around the boundaries have undefined curvature, calling for use of certain constraints in the computation, without which the boundaries would have corrupted the caricature.

FIGS. 4 and 5, namely, demonstrate how regions with high Gaussian curvature scale up, whereas regions with small Gaussian curvature shrink. Regions with approximately constant Gaussian Curvature were scaled uniformly. These caricatures were created using the geometric driven approach. That is, without a reference model.

We tested our method on faces and shapes with boundaries, see FIGS. 5, 6, and 7. The resulting caricatures demonstrate an exaggeration of unique features and expressions of the subject's face in a natural comical fashion. The nose was stretched while flat regions like the cheeks, shrank, suggesting the applicability of our method to face recognition.

In FIG. 6, from left to right: an original facial 3D model having a texture; a caricature with an exaggeration factor γ=0.5; and a caricature with an exaggeration factor γ=0.75.

In FIG. 7, from left to right: another original facial 3D model; a caricature with an exaggeration factor γ=0.5; and a caricature with an exaggeration factor γ=0.75.

As with most mesh deformation techniques, the present method may introduce self-intersections, particulately for large exaggeration factors. However, since we solve a Poisson equation, this intersection occurs only globally. Using post-processing self-collision handling methods, such as the one of D. Harmon, D. Panozzo, O. Sorkine, D. Zorin, “Interference-aware geometric modeling”, in ACM Transactions on Graphics (TOG), Vol. 30, ACM, 2011, p. 137, one can remove these intersections.

3.2 Run-Time Analysis

The algorithm was implemented both in Matlab and in C++. For the C++ version, we used the open source version of Eigen Library (G. Guennebaud, B. Jacob, et al., “Eigen v3”, http://eigen.tuxfamily.org, 2010, last viewed Apr. 5, 2015). All the experiments were executed on a 2.6 GHz Intel Core i7 machine with 8 Gigabytes (GB) of Random-Access Memory (RAM). Run-times for various meshes from the TOSCA dataset (A. Bronstein, M. Bronstein, R. Kimmel, “Project TOSCA”, 2010, http://tosca.cs.technion.ac.il, last viewed Apr. 5, 2015) are shown in Table 1. The preprocessing stage includes calculating the laplacian matrix, the Gaussian curvature, the gradient of the original shape, and matrices for calculating the discrete divergence efficiently. The core runtime includes taking the power of the curvature in each triangle, multiplying the gradients fields by |K|^(γ), and solving the resulting sparse linear system.

TABLE 1 Pre-processing and core run-times (in milliseconds) of the present algorithm for exaggerating features of shapes from TOSCA dataset. Model No. of No. of Pre-processing Core name vertices faces runtime runtime Wolf 4344 8684 349.5 7.4 Centaur 15768 31532 1538.2 29.3 Dog 25290 50528 2944.7 47.0 Human 52565 105032 7245.6 105.8

3.3 Robustness to Almost Isometric Deformation Analysis

To evaluate the robustness of the present method to almost isometric deformations of shapes, we applied it to shapes in two different poses. We calculated the histograms of normalized distortion of geodesic distances between corresponding points on the two shapes, before and after exaggeration was applied (see FIGS. 8 and 9). The results demonstrate that caricaturizing two almost isometric shapes indeed produces two almost isometric caricatures.

FIG. 8 shows a histogram of relative geodesic distortion (in %) between caricatures of a Centaur. The X axis denotes the amount of point pairs, and the Y axis denotes the percent change in geodesic distance between these two points. From bottom to top: Almost isometric centaurs; centaurs with exaggeration factor γ=0.2; and centaurs with exaggeration factor γ=0.3.

FIG. 9 shows a histogram of relative geodesic distortion (in %) between caricatures of a cat. The X axis denotes the amount of point pairs, and the Y axis denotes the percent change in geodesic distance between these two points. From bottom to top: Almost isometric cats; cats with exaggeration factor γ=0.2; cats with exaggeration factor γ=0.3; and a map which is as isometric as possible between a cat and a dog.

To demonstrate the isometry preservation property, we exaggerated a centaur shape in two different poses. In FIG. 10, top row, from left to right: an original 3D model of the centaur in a first pose; a caricature of the first pose with an exaggeration factor γ=0.25; and a caricature of the first pose with an exaggeration factor γ=0.5. Bottom row, left to right: an original 3D model of the centaur in a second pose; a caricature of the second pose with an exaggeration factor γ=0.25; and a caricature of the second pose with an exaggeration factor γ=0.5.

FIGS. 11, 12 and 13 show examples of exaggerated shapes in various poses for different values of γ.

In FIG. 11, the top row shows original 3D models of a dog in four different poses. The bottom row shows caricatures of these four different poses with an exaggeration factor γ=0.3.

In FIG. 12, the top row shows original 3D models of a gymnast in four different poses. The middle row shows caricatures of these four different poses with an exaggeration factor γ=0.3. The bottom row shows caricatures of these four different poses with an exaggeration factor γ=0.4.

In FIG. 13, the top row shows original 3D models of a martial artist in four different poses. The middle row shows caricatures of these four different poses with an exaggeration factor γ=0.3. The bottom row shows caricatures of these four different poses with an exaggeration factor γ=0.4.

The Intrinsic Surface Caricaturization algorithm was successfully harnessed for the present audio-based caricaturization method. Artificial motion for the video sequence was generated by linearly interpolating between different caricaturized poses of the same object. This interpolation was then projected onto the space of isometrics by solving a Poisson equation, as described in M. Botsch, R. Sumner, M. Pauly, M. Gross, “Deformation transfer for detail-preserving surface editing”, in Vision, Modeling & Visualization, Citeseer, 2006, 357-364, and in R. W. Sumner, J. Popovic, “Deformation transfer for triangle meshes”, in ACM Transactions on Graphics (TOG), Vol. 23, ACM, 2004, 399-405.

3.4 Shape Classification by Exaggeration

A prominent task in shape analysis is the ability of a computer to classify shapes. To facilitate this task, we have exaggerated the features of the query shapes as a preprocessing step before classification. This method is justified by an experiment which showed that humans perform better in memorizing and recognizing caricaturized faces than regular ones. That experiment is described in R. Kurzweil, “How to create a mind: The secret of human thought revealed”, Penguin, 2012. Kurzweil mentions the ability of the brain to recognize the subject behind a caricature. That is, the human pattern recognition ability is robust to some classes of exaggerations.

We have successfully demonstrated the power of caricaturisation in facilitating recognition of 3D shapes that have undergone various almost-isometric deformations. First, we manually sampled points on quadruped mammals in various poses. These points included, for example, the tip of the tail, the tip of the nose, etc. of the mammal We then normalized the shapes and calculated the geodesic distances between every two points, forming a matrix sized approximately 30×30 for each shape. Each matrix was flattened into an approximately 900-long vector. We then calculated the L2 distance between every two vectors, and placed the results in a matrix D. In order to represent each shape as a point in the plane, we calculated the Multidimentional Scaling of D to the plane. That is, we find the configuration of points in the plane between which the distances are as close as possible (in least squares sense) to the L2 distances between the representations of the shapes in the higher dimensional (900-long) space. We repeated the process after exaggerating the shapes.

The embeddings of the shapes to the plane are presented in FIGS. 14A-B and 15A-B. The exaggeration is shown to separate the mammal classes, making the classification task simpler.

FIG. 14A shows a graph of the results of the Multidimentional Scaling for uncaricaturized versions of a dog, a cat and a wolf Notably, a dog area 100 in the graph and a wolf area in the graph 104 greatly overlap, whereas a cat area 102 is relatively distinct from the previous two. This means that, in these uncaricaturized 3D models, distingushing between models using the Multidimentional Scaling of the geodesic distance matrices would yield unsatisfactory results.

FIG. 14B shows a similar graph, this time for caricaturized versions (γ=0.3) of the dog, cat and wolf. As clearly seen in this figure, a dog area 200, a cat area 202 and a wolf area 204 are easily distinguishable using simple cartesian analysis of this graph. Here, the Multidimentional Scaling of the geodesic distance matrices has proven efficient in allowing for computerized classification of 3D shapes.

FIG. 15A is similar to the graph of FIG. 14A, but pertains to an added quadruped mammal—a horse. Here, again, distingushing between the shapes is not feasible—especially between the dog, wolf and horse.

FIG. 15B shows the same dog, cat, wolf and horse—this time after caricaturization (γ=0.3). The Multidimentional Scaling of the geodesic distance matrices has proven efficient also here.

4.5 Application: Exaggeration Volume

The applicability of the Intrinsic Surface Caricaturization algorithm to the present audio-based caricaturization method was successfully tested, and a number of audiovisual clips were synthesized. We first constructed a look-up table of caricatures for equally sampled values of the exaggeration factor γε[0,1]. For each video frame, we calculated the mean absolute value volume of the audio within the corresponding time-frame of the audio sequence. Then, we found its the closest exaggeration factor γ in the look-up table and set the video frame to its corresponding rendered caricature. The resulting audiovisual clip showed a 3D shape caricaturized in synchronization with the played song.

The present invention may be a system, a method, and/or a computer program product. The computer program product may include a computer readable storage medium (or media) having computer readable program instructions thereon for causing a processor to carry out aspects of the present invention.

The computer readable storage medium can be a tangible device that can retain and store instructions for use by an instruction execution device. The computer readable storage medium may be, for example, but is not limited to, an electronic storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the foregoing. A non-exhaustive list of more specific examples of the computer readable storage medium includes the following: a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), a static random access memory (SRAM), a portable compact disc read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, a mechanically encoded device, or any suitable combination of the foregoing. A computer readable storage medium, as used herein, is not to be construed as being transitory signals per se, such as radio waves or other freely propagating electromagnetic waves, electromagnetic waves propagating through a waveguide or other transmission media (e.g., light pulses passing through a fiber-optic cable), or electrical signals transmitted through a wire.

Computer readable program instructions described herein can be downloaded to respective computing/processing devices from a computer readable storage medium or to an external computer or external storage device via a network, for example, the Internet, a local area network, a wide area network and/or a wireless network. The network may comprise copper transmission cables, optical transmission fibers, wireless transmission, routers, firewalls, switches, gateway computers and/or edge servers. A network adapter card or network interface in each computing/processing device receives computer readable program instructions from the network and forwards the computer readable program instructions for storage in a computer readable storage medium within the respective computing/processing device.

Computer readable program instructions for carrying out operations of the present invention may be assembler instructions, instruction-set-architecture (ISA) instructions, machine instructions, machine dependent instructions, microcode, firmware instructions, state-setting data, or either source code or object code written in any combination of one or more programming languages, including an object oriented programming language such as Java, Smalltalk, C++ or the like, and conventional procedural programming languages, such as the “C” programming language or similar programming languages. The computer readable program instructions may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the latter scenario, the remote computer may be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet Service Provider). In some embodiments, electronic circuitry including, for example, programmable logic circuitry, field-programmable gate arrays (FPGA), or programmable logic arrays (PLA) may execute the computer readable program instructions by utilizing state information of the computer readable program instructions to personalize the electronic circuitry, in order to perform aspects of the present invention.

Aspects of the present invention are described herein with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and/or block diagrams, and combinations of blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer readable program instructions.

These computer readable program instructions may be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks. These computer readable program instructions may also be stored in a computer readable storage medium that can direct a computer, a programmable data processing apparatus, and/or other devices to function in a particular manner, such that the computer readable storage medium having instructions stored therein comprises an article of manufacture including instructions which implement aspects of the function/act specified in the flowchart and/or block diagram block or blocks.

The computer readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable apparatus or other device to produce a computer implemented process, such that the instructions which execute on the computer, other programmable apparatus, or other device implement the functions/acts specified in the flowchart and/or block diagram block or blocks.

The descriptions of the various embodiments of the present invention have been presented for purposes of illustration, but are not intended to be exhaustive or limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments. The terminology used herein was chosen to best explain the principles of the embodiments, the practical application or technical improvement over technologies found in the marketplace, or to enable others of ordinary skill in the art to understand the embodiments disclosed herein. 

We claim:
 1. A method comprising using at least one hardware processor for: receiving a three-dimensional model of an object, wherein the three-dimensional model is embodied as a digital file that comprises a representation of the object; receiving an audio sequence embodied as a digital file that comprises a musical composition; generating a video frame sequence, wherein the generating comprises computing a caricature of the object by applying a computerized caricaturization algorithm to the three-dimensional model, wherein the computing comprises: scaling gradient fields of surface coordinates of the three-dimensional model by a function of a Gaussian curvature of the surface, and finding a regular surface whose gradient fields fit the scaled gradient fields, wherein (a) the computing is with a different exaggeration factor for each of multiple ones of the video frames, and (b) the different exaggeration factor is based on one or more parameters of the musical composition of the audio sequence; and synthesizing the audio sequence and the video frame sequence into an audiovisual clip.
 2. The method according to claim 1, further comprising using said at least one hardware processor for determining the one or more parameters for each of multiple periods of the audio sequence, wherein the one or more parameters are selected from the group consisting of: amplitude, frequency and tempo.
 3. The method according to claim 1, wherein the generating further comprises altering a view angle of the caricature along the video frame sequence.
 4. The method according to claim 1, wherein the exaggeration factor is applied uniformly, to the entirety of the three-dimensional model.
 5. The method according to claim 1, wherein the exaggeration factor is applied non-uniformly, only to one or more portions of the three-dimensional model, which portions amount to less than the entirety of the three-dimensional model.
 6. The method according to claim 1, wherein the computing of the caricature of the object further comprises: constructing a look-up table comprised of (a) different visualizations of the caricature, each computed with one of the different exaggeration factors, and (b) the exaggeration factor for each of the different visualizations; and using each caricature visualization from the look-up table when the exaggeration factor of that caricature visualization is determined to be suitable for the one or more parameters of the audio sequence.
 7. The method according to claim 1, wherein the computing of the caricature of the object further comprises amplifying the scaling according to local discrepancies between the object and a reference object, wherein the reference object is not a scaled down version of the object.
 8. The method according to claim 1, wherein the computing of the caricature of the object further comprises amplifying the scaling according to local discrepancies between the object and a scaled down version of the object.
 9. A computer program product comprising a non-transitory computer-readable storage medium having program code embodied thereon, the program code executable by at least one hardware processor for: receiving a three-dimensional model of an object, wherein the three-dimensional model is embodied as a digital file that comprises a representation of the object; receiving an audio sequence embodied as a digital file that comprises a musical composition; generating a video frame sequence, wherein the generating comprises computing a caricature of the object by applying a computerized caricaturization algorithm to the three-dimensional model, wherein the computing comprises: scaling gradient fields of surface coordinates of the three-dimensional model by a function of a Gaussian curvature of the surface, and finding a regular surface whose gradient fields fit the scaled gradient fields, wherein (a) the computing is with a different exaggeration factor for each of multiple ones of the video frames, and (b) the different exaggeration factor is based on one or more parameters of the musical composition of the audio sequence; and synthesizing the audio sequence and the video frame sequence into an audiovisual clip.
 10. The computer program product according to claim 9, wherein: the program code is further executable by said at least one hardware processor for determining the one or more parameters for each of multiple periods of the audio sequence; and the one or more parameters are selected from the group consisting of: amplitude, frequency and tempo.
 11. The computer program product according to claim 9, wherein the generating further comprises altering a view angle of the caricature along the video frame sequence.
 12. The computer program product according to claim 9, wherein the exaggeration factor is applied uniformly, to the entirety of the three-dimensional model.
 13. The computer program product according to claim 9, wherein the exaggeration factor is applied non-uniformly, only to one or more portions of the three-dimensional model, which portions amount to less than the entirety of the three-dimensional model.
 14. The computer program product according to claim 9, wherein the computing of the caricature of the object further comprises: constructing a look-up table comprised of (a) different visualizations of the caricature, each computed with one of the different exaggeration factors, and (b) the exaggeration factor for each of the different visualizations; and using each caricature visualization from the look-up table when the exaggeration factor of that caricature visualization is determined to be suitable for the one or more parameters of the audio sequence.
 15. The computer program product according to claim 9, wherein the computing of the caricature of the object further comprises amplifying the scaling according to local discrepancies between the object and a reference object, wherein the reference object is not a scaled down version of the object.
 16. The computer program product according to claim 9, wherein the computing of the caricature of the object further comprises amplifying the scaling according to local discrepancies between the object and a scaled down version of the object.
 17. A system comprising: (a) a non-transitory computer-readable storage medium having program code embodied thereon, the program code comprising instructions for: receiving a three-dimensional model of an object, wherein the three-dimensional model is embodied as a digital file that comprises a representation of the object, receiving an audio sequence embodied as a digital file that comprises a musical composition, generating a video frame sequence, wherein the generating comprises computing a caricature of the object by applying a computerized caricaturization algorithm to the three-dimensional model, wherein the computing comprises: scaling gradient fields of surface coordinates of the three-dimensional model by a function of a Gaussian curvature of the surface, and finding a regular surface whose gradient fields fit the scaled gradient fields, wherein (i) the computing is with a different exaggeration factor for each of multiple ones of the video frames, and (ii) the different exaggeration factor is based on one or more parameters of the musical composition of the audio sequence, and synthesizing the audio sequence and the video frame sequence into an audiovisual clip; and (b) at least one hardware processor configured to execute the instructions.
 18. The system according to claim 17, wherein: the program code is further executable by said at least one hardware processor for determining the one or more parameters for each of multiple periods of the audio sequence; and the one or more parameters are selected from the group consisting of: amplitude, frequency and tempo.
 19. The system according to claim 17, wherein the generating further comprises altering a view angle of the caricature along the video frame sequence.
 20. The system according to claim 17, wherein the exaggeration factor is applied uniformly, to the entirety of the three-dimensional model.
 21. The system according to claim 17, wherein the exaggeration factor is applied non-uniformly, only to one or more portions of the three-dimensional model, which portions amount to less than the entirety of the three-dimensional model.
 22. The system according to claim 17, wherein the computing of the caricature of the object further comprises: constructing a look-up table comprised of (a) different visualizations of the caricature, each computed with one of the different exaggeration factors, and (b) the exaggeration factor for each of the different visualizations; and using each caricature visualization from the look-up table when the exaggeration factor of that caricature visualization is determined to be suitable for the one or more parameters of the audio sequence.
 23. The system according to claim 17, wherein the computing of the caricature of the object further comprises amplifying the scaling according to local discrepancies between the object and a reference object, wherein the reference object is not a scaled down version of the object.
 24. The system according to claim 17, wherein the computing of the caricature of the object further comprises amplifying the scaling according to local discrepancies between the object and a scaled down version of the object. 